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ABSTRACT 

Coalescing binary black holes experience an impulsive kick due to anisotropic emission of gravitational 
waves. We discuss the dynamical consequences of the recoil accompanying massive black hole mergers. 
Recoil velocities are sufficient to eject most coalescing black holes from dwarf galaxies and globular 
clusters, which may explain the apparent absence of massive black holes in these systems. Ejection 
from giant elliptical galaxies would be rare, but coalescing black holes are displaced from the center and 
fall back on a time scale of order the half-mass crossing time. Displacement of the black holes transfers 
energy to the stars in the nucleus and can convert a steep density cusp into a core. Radiation recoil 
calls into question models that grow supermassive black holes from hierarchical mergers of stellar-mass 
precursors. 

Subject headings: black hole physics — gravitation — gravitational waves — galaxies: nuclei 



1. KICK AMPLITUDE 

In a companion paper l|Favata. Hughes fc Hol3l200l 
hereafter Paper I) , the amplitude of the recoil velocity ex- 
perienced by a binary black hole (BH) due to anisotropic 
emission of gravitational radiation during coalescence is 
computed . Here we explore some o f the consequences of 
the kicks (Rcdmount & Rces 198|j): the probability that 
BHs are ejected from galaxies and the implications for 
BH growth; the time scale for a kicked BH to return to 
the center of a galaxy; the effect of displacement on nu- 
clear structure; and other observational signatures of the 
kicks. Unless otherwise indicated, notation is the same 
as in Paper I. 

For inspiral from a circular orbit, the kick velocity is a 
function of the binary mass ratio q = ?tii/to2 < 1, the BH 
spins 5i and ct2, and the initial angle l between the spin of 
the larger BH and the orbital angular momentum of the 
binary. Following Paper I, the spin of the smaller BH is 
ignored. Although Paper I only considers the cases t = 
and l — 180, the recoil for arbitrary inclination is likely 
to be bounded between these extreme values. Also, the 
detailed inclination dependence is unimportant in com- 
parison with the large uncertainty already present in the 
contribution to the recoil from the final plunge and coa- 
lescence. We will therefore assume that the restriction to 
equatorial- prograde/retrograde orbits (5,2 = [—1, 1]) con- 
sidered in Paper I encompasses the characteristic range 
of recoil velocities. 

Figure 2b of Paper I shows upper- and lower-limit esti- 
mates of the recoil velocity as a function of the effective 
spin parameter 5, for reduced mass ratio i] = fJ./M = 
g/(l + q) 2 = 0.1. The upper limit for r\ = 0.1 is well fit 
in the range —0.9 < a < 0.8 by the following fifth-order 
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polynomial: 

Kipper = 465 km s _1 ^-(l - 0.2815 - 0.03615 2 

/max 

-0.3465 3 - 0.3745 4 - 0.1845 5 ). (1) 

Fitchett's (1983) scaling function /(g)// m ax, with /(g) = 
g 2 (l - q)/(l + g) 5 , equals 0.433 for r) = 0.1. The lower 
limit curve of Paper I is well fit by 

U lowcr = 54.4 km s" 1 ^-(l + 1.225+ 1.04a 2 

/max 

+ 0.9775 3 - 0.2015 4 - 0.4345 5 ). (2) 

We convert these expressions into estimates of the 
bounds on Vkick as follows. First, as discussed in Paper 
I, there is an ambiguity in how one translates the physi- 
cal spin parameter 52 of the larger hole into the effective 
spi n parameter 5 of equations Q and (J2J) . Here we adopt 
the iDamourl l)200H) relation a = (1 + 3g/4)(l + g)~ 2 5 2 . 
Second, Fitchett's scaling function assumes that both 
bodies are non-spinning, and vanishes when q = 1. In 
fact, when 5^0, significant recoil would occur even 
for q = 1 due to spin-orbit coupling. We can guess the 
approximate form of a new s caling functio n by examin- 
ing the spin-orbit corrections l)Kidderlll995l) to Fitchett's 
recoil formula. For equatorial orbits, equation (4) of Pa- 
per I suggests that /(g) should be multiplied by the fac- 
tor |1 + (7/29)22/(1 - g)|/|l + (7/29)5 2 /(l - g')|, where 
g' = 0.127 is the value used in defining V uppor and Ui OW er 
in equations Jpl and (J2J). 

Figure ^ plots upper and lower limits to Vkick as func- 
tions of 52 and q. The average over 52 of the up- 
per limit estimates are ~ (138,444,154) km s _1 for 
q = (0.1,0.4,0.8); Figure H suggests a weak dependence 
on a2- Lower limit estimates are more strongly spin- 
dependent; the averages over 0,2 are ~ (21.1,63.6,24.9) 
km s _1 for the same values of q. For moderately large 
spins (52 > 0.8) and prograde capture, the lower limit 
estimates exceed 100 km s _1 for 0.2 < q < 0.6. In what 
follows, we will assume that ~ 500 km s _1 is an absolute 
upper limit to Vkick- 

2. ESCAPE 
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When Vkick > V csc = ^/2<p(r = 0), with </>(r) the grav- 
itational potential of the system (galaxy, dark matter 
halo) hosting the BH, the BH has enough kinetic energy 
to escape. Figure|21shows central escape velocities in four 
types of stellar system that could contain merging BHs: 
giant elliptical galaxies (E), dwarf ellipticals (dE), dwarf 
spheroidals (dSph) and globular clusters (GC). We fit the 
trend log (y esc /l km s -1 ) = A — /3My separately for each 
class of object. dEs and GCs each separately establish a 
relation L ~ V^ Lj for GCs, this is comp atible with the 
relation found bv lDiorgovski et ahl (^997). The E sample 
is consistent with the Faber- Jackson (1976) relation. 
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Fig. 1. — Upper limit (top) and lower limit (bottom) es- 
timates of Vkick as functions of mass ratio q and spin of the 
larger black hole 5,2- Units are km s — 1 . Values of a 2 and q cor- 
responding to a > 0.8 lie in the region to the right of the dotted 
line. Since equations Q and J2J are not valid for a > 0.8, a 
was replaced by 0.8 in this region. 

The solid line in Figure shows escape velocities from 
the dark matter (DM) halos associated with the lumi- 
nous stellar systems. To relate halo properties to galaxy 
luminosities, we use the conditional luminosity function 
$(L\M)dL from the concordance ACDM model Ml of 
lYang. Mo k van den Boscbl l)2003|) . This function gives 
the number of galaxies in the luminosity range L ± dL/2 
that reside inside a halo with virial mass Af V j r . The av- 
erage luminosity L\ of the brightest ("central") galaxy 
in the halo of mass M v ; r is implicitly given by the con- 
dition J L §{L\M v i r )dL = 1. Inverting this we ob- 
tain M v i r (Iri) and relate this mass to the escape ve- 
locity via V^ sc — 2cg(c)GM v i T /R v i T where i? v ir is the 
virial radius of the halo, c is the concentration of a 
halo obeying the Navarro, Frenk k White (1996; here- 
after NFW) profile, and q(c ) = [ln(l + c) - c/(l + c)}" 1 
(e.g.. lLokas k Mamonl l200l). Both R y-, r and c are func- 
tions of Myj r and the red shift z (e.g.. iBrvan k Normanl 
119981 iBullock et all 1200 lj) . At z = the average es- 
cape velocity is given by V csc = 239 km s~ (mn/h) 1 ^ 2 , 
where M V - 1T = lO^mnMg) and h is the Hubble parame- 
ter, set to 0.7 in Figure|21 The upturn in escape velocity 
for galaxies brighter than My ~ —20 is a consequence 



of the increase in the occupation number of their host 
halos. The dashed line in Figure El shows the escape ve- 
locity from the combined luminous 4- DM potential for 
the E galaxies, using the scaling relation derived above 
to describe the luminous component. 

Figure suggests that the consequences of the kicks 
are strikingly different for the different classes of stel- 
lar system that might host BHs. Escape velocities 
from E galaxies are dominated by the stellar contri- 
bution to the potential; in the sample of Faber et 
al. (1997), V csc > 450 km s^ 1 even without accounting 
for DM. This exceeds even the upper limits in Figure 
Hence, the kicks should almost never unbind BHs 
from E galaxies. The tight corr elations observed between 
BH mass and bulge luminosity jMcClure k Dunlorl2002t 
lErwin. Graham k. Caonl I2003D and velocity disp ersion 
(jFerrarese k MerrittJ 120001: iGebhardt et alJlSooqi could 
probably not be maintained if escape occurred with any 
significant frequency from luminous galaxies. 

The existence of DM significantly affects the escape 
probability from dE and dSph galaxies, implying kicks 
of ~ 300 km s _1 and ~ 100 km s _1 respectively for es- 
cape. In the absence of DM, these numbers would be 
~ 100 km s" 1 and ~ 20 km s~ 1 respectively. Hence, 
kicks of order 200 km s _1 would unbind BHs from dSph 
galaxies whether or not they contain DM, while dE galax- 
ies could retain their BHs if they are surrounded by DM 
halos. 
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Fig. 2. — Central escape velocities in km s" 1 in four 
types of stellar syste m that coul d harb or merging BHs. E 
galaxy data are from Faber et al. (1997), with separate sym- 
bols for core (□) and pow er-law (A) galaxies. dE data a re from 
Binuecli & r, Jerierj I1998A . with mass-to-light ratios from Mateo 
1 1998). GC and dSph data are from the tabulation of Wcbbink 
( 1996). Solid line is the mean escape velocity from the DM halos 
associated with the luminous matter. Dashed line is the escape 
velocity from the combined luminous + mean DM potentials 
for E galaxies. 

Evidence for intermediate-mass black holes at the cen- 
ters of galaxies fainter t han My sw —19 is sketchy 
(e.g., Ivan der Marell l2003|) . although there is indirect 
non-dynamical) eviden ce for BHs in faint Seyfert bulges 
Filippenko k Holl20 03). We note that the dense nuclei 
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associated with BHs in galaxies like M32 (My w —19) 
become progressively less frequent at magnitudes fainter 
than My ps — 16 and d isappear entirely below My w — 12 
ijvan den Berghl Il986j) . If the dense nuc lei are associ- 
ated with nuclear BHs (e.g.. lPeeblesHl972Tl . their absence 
could signal loss of the BHs via ejection. It is intriguing 
that these nuclei are someti mes observed to be displaced 
far fr om the galaxy center (|Binggeli. Barazza fe Jerienl 
2000) . Figures and [21 imply that even kicks at the 
lower limits of Paper I would almost always unbind BHs 
from GCs. 

3. EJECTION IN HIERARCHICAL MERGING SCENARIOS 

The kicks have serious implications for models in which 
massive BHs grow from mergers of less massive seeds. 
In some of these models, the precursors are stellar- 
or intermediate-mass black holes produced in the col- 
lapse of the first stars (Population III) and the merg- 
ing c ommenced in minihalos at redshifts as large as 
~ 20 ilMadau fe Rees»200lHyolonteri. Haardt. fe Mad"aul 
120031 llslam. Tavlor fe Silk! 12003(1 . We evaluate the 
plausibility of such models in light of the estimates 
of Vkick derived in Paper I. Kicks from gravitational 
wave emission may compete with high-velocity recoils 
( Sasl aw. Valtonen. fc Aarsethl Il974|) from (Newtonian) 
three-body interactions. While the Newtonian recoil oc- 
curs only when three BHs are present, which is contin- 
gent on the galaxy merger rate and the BH binary orbital 
decay rate, radiation recoil is present whenever BHs co- 
alesce. 
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Fig. 3. — The maximum redshift z c j Gct at which (a) DM halos 
only, and (b) DM halos and the central galaxies combined, can 
confine BHs as a function of the z = BH mass, for five values of 
the kick velocity. The depth of the galactic contribution to the 
potential was calculated by identifying the velocity dispersion 
of the stellar sp heroid with the circular velocity of the halo 
IFerraresdl200ai . 

The confining effect of DM halos in a hierarchical uni- 
verse was smaller at higher redshifts when the average 
halo mass was smaller. We estimate the maximum red- 
shift at which DM h alos can confine the progenitors of 
the present-day BHs. iFerraresel (|2002|) derived a relation 
of the present-day BH mass M^(z = ) to the mass of 
the ho st halo M vir (z = 0). We use the IWvithe fc Loebl 
(2002) form of the relation (their equation 11) to obtain 
the host halo mass and extrapolate t he mass back in red- 
shift via the accretion history mode l l)Bullock et all20oHl 
calibrated bv lWechsler et alJ <|2002f) on a set of numerical 
simulations of DM clustering in a ACDM universe. The 
accretion trajectory M v iy-(z) oc e~ az , where a is itself a 
function of the halo mass at z = 0, can be interpreted as 
the mass of the most massive, and thus the most easily 
confining parent halo at redshift z. We can then calculate 



the escape velocity V esc (z) of the most massive progen- 
itor halo as a function of redshift. Finally, we solve for 
z ojcct such that Vkick = K> sc Reject); this is the maximum 
redshift at which the progenitors of the present-day BHs 
could have started merging. We also modelled the ef- 
fect on Zcjoct of including the potential due to a stellar 
component, idealized as an isothermal sphere with core 
radius 77, = 2GM^/a 2 and outer cutoff 10 3 r^,; the ID 
stellar velocity dispersion a was related to halo circular 
velocity as in Ferrarese (2002). 

The results for five representative choices of Vkick are 
shown in Figure For Vkick ~ 150 km s~\ we find 
Reject < 11(14) over the entire range of Mbh\ the latter 
value is from the models that include a stellar compo- 
nent. For Vkick ~ 300 km s^ 1 , the assembly of a 10 8 A/ Q 
BH must have started at z < 8(10). Models that grow 
supcrmassive BHs from mergers of seeds of much lower 
mass at redshifts z > 10 are thus disfavored due to the 
difficulty of retaining the kicked BHs. The effects of the 
kicks could be mitigated if early growth of the BHs was 
dominated by accretion, or if BH mergers were delayed 
until their halos had grown much more massive. 

4. FALLBACK TIMES 

A BH that has been kicked from the center of a stel- 
lar system with a velocity less than V csc falls back and 
its orbit decays via dynamical friction against the stars 
and gas. We define the fallback time Ti n f a n as the time 
required for a BH to return to a zero-velocity state after 
being ejected. The velocity with which the BH is ejected 
from the site of the merger is Kject = (Mbh/M e g) Vkick < 
Vkick; here M e & = M hh + M hound with M bound the 
mass in stars that remain bound to the BH after it is 
kicked. For recoil in a singular isothermal sphere nu- 
cleus p oc r- 2 , Mcff/Afbh « (1.9,1.5,1.05,1.00) when 
Vkick/c = (0.5, 1, 2, 3) where a is the ID stellar velocity 
dispersion; M bo und/M b h oc (Vi ick /cr) -4 for Vkick > cr. 
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Fig. 4. — Effect on the nuclear density profile of black hole 
ejection. The initial galaxy model (black line) has a p ~ r~ 1 
density cusp, (a) Impulsive removal of the black hole. Tick 
marks show the radius of the black hole's sphere of influence 
before ejection. A core forms with radius ~ 2r^. (b) Ejection at 
velocities less than V CB c- The black hole has mass 0.003Af ga i; the 
galaxy is initially spherical and the black hole's orbit remains 
nearly radial as it decays via dynamical friction. The arrow in 
panel marks r%. 

We evaluated Ti n f a n for BHs kicked from the centers 
of Dehnen (1993)-law galaxies for which the central den- 
sity obeys p oc r -7 . Bright E galaxies have < 7 < 1 
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ijGebhardt et al.lH996j) . and cusps steeper than this are 
likely to be softened by the binary BH prior to coales- 
cence llMilosavlievic fc MerrittJl2~001fl and by the ejection 
itself (§5). Given values for M c g and V^j ec t, the fall- 
back time in a spherical galaxy is gi ven by the orbit- 
avera ged dynamical friction equation ijCohn fc Kulsrudl 
fl978h . For V e j cct /V csc < 0.6, infall times were found to be 

well approximated by T infall rj T 1/2 (V r C j OCt /V; s c) 2 - 5(1+7) 
for M eff = 0.001Mg a i, where the period Ti/ 2 of a cir- 
cular orbit at the galaxy's half-light radius is given in 
terms of the galaxy's visual luminosity by T\/2 ~ 2 x 

10 s yr^y/lO 11 ^) 1 / 2 llValluri fc Merritd Il99^ . Thus, 
return of a BH to a stationary state requires of order a 
few times 10 s yr or less over a wide range of cusp slopes 
and galaxy luminosities for V e ^ ec t V csc /2. As indicated 
in Figure [21 this is the likely situation in the bright E 
galaxies. Infall times are especially short for 7 > 1, since 
the BH experiences a strong impulsive frictional force 
as it passes repeatedly through the dense center. When 
Reject ;$ C) the BH never moves far from its central po- 
sition and it carries much of the nucleus with it. We 
carried out A-body simulations of this regime and found 
that return to zero velocity occurs in roughly one orbital 
period when V^ct ^ In fainter dE and dSph galax- 
ies, ejection would more often occur near V osc and infall 
times could be arbitrarily long, determined primarily by 
the mass distribution at large radii. 

In a nonspherical galaxy, an ejected BH does not pass 
precisely through the dense center on each return, de- 
laying the infall. To test the effect of non-spherical 
geometries on the infall time, we carried out experi- 
me nts in the triaxial generali zations of the Dehnen mod- 
els l)Merritt fc Fridmanl 1996). Results were found to de- 
pend only weakly on the axis ratios of the models. Decay 
times in the triaxial geometry exhibit a spread in values 
depending on the initial launch angle, bounded from be- 
low by the decay time along the short axis. We found a 
mean at every V e ject/V e sc that is ~ 3 — 5 times greater 
than in a spherical galaxy with the same cusp slope. 

5. OBSERVABLE CONSEQUENCES OF THE 
DISPLACEMENT 

Displacement of the BH also transfers energy to the 
nucleus and lowers its density within a region of size ~ 
Th, the radius of the BH's sphere of influence (defined 
here as the radius of a sphere containing a mass in stars 
equal to twice that of the BH). The simplest case to 
consider is I4j ec t ^ 14 sc ; the BH and its entrained mass 
depart the nucleus on a time scale that is of order the 



crossing time at r/j or less and do not return. The effect 
on the nucleus can be approximated by constructing a 
steady-state model of a galaxy containing a central point 
mass, then removing the point mass instantaneously and 
allowing the remaining particles to relax to a new steady 
state. Figure 0^ shows the results for three values of 
M e s/M ga \. Initial conditions consisted of 10 6 particles 
representing stars in a 7 = 1 Dehnen model. We find 
that a core of roughly constant density forms within a 
radius of ~ 2r/j. Setting 7 = 2 (not shown) results in 
a core of size ~ rv Figure^ shows the change in the 
nuclear density profile for simulations with V^ct < Vesc- 
Significant changes in the central density require V^j CC t > 
0.25V CSC . We conclude that the recoil could affect the 
observable structure of nuclei, sinc e radii of ~ Ith are 
resolv ed in many nearby galaxies l|Merritt fc Ferraresel 

2001) . 

The "mass deficits" seen at the c enters 
of bright galaxies llMilosavlieyic et alJ l2002t 
iRavindranath. Ho fc Filippenkol 120021) may be due 
to the combined effects of slingshot ejection and BH 
displacement, although we note that the large cores 
observed in some bright galaxies c ould probably not 
be p roduced by either mechanism ( Mil osavlievic et all 

2002) . 

The X -shaped radio sources associated with giant E 
galaxies (|Leahv fc Parmall 19921) are plausible sites of re- 
cent BH coalesceiicelMe^^^^Ete3llQn2)- Displace- 
ment of the merged BHs from the galaxy center prior to 
ignition of the "active" lobes would imply a distortion of 
the X-morphology, in the sense that the "wings" (the in- 
active lobes) would be non-collinear near the center of the 
X. Such di stortions are in fact a common feature of the 
X-sources ijGopal- Krishna. Biermann fc Wiitall2003[) . al- 
though the linear scale of the distorti ons in some of the X- 
sources (e.g. ~ 10 kpc in NGC 236: iMurgia et alJl2001l) 
suggests that orbital motion o f the merging galaxi es may 
be a more likely explanation jBalcells et al.lll9 95^ . 
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